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Abstract 

We study the initial and final state interaction effects in the transverse momentum dependent 
parton distributions in the small- £ saturation region. In particular, we discuss the quark dis- 
tributions in the semi-inclusive deep inelastic scattering, Drell-Yan lepton pair production and 
dijet-correlation processes in pA collisions. We calculate the quark distributions in the scalar-QED 
model and then extend to the color glass condensate formalism in QCD. The quark distributions 
are found universal between the DIS and Drell-Yan processes. On the other hand, the quark dis- 
tribution from the qq' —> qq' channel contribution to the dijet-correlation process is not universal. 
However, we find that it can be related to the quark distribution in DIS process by a convolution 
with the normalized unintegrated gluon distribution in the color glass condensate formalism in the 
large N c limit. 
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I. INTRODUCTION 



Partonic internal structure of nucleon and nucleus have attracted many theoretical and 
experimental investigations in the past and are still in the frontier of the subatomic physics 
research. These studies aim at providing us accurate description of the hadronic structure in 
terms of fundamental degree of freedom in Quantum Chromodynamics (QCD), and mean- 
while presenting an important path to discover the new physics beyond the Standard Model, 
which are currently undertaken at various high energy experiments, such as the Fermilab 
Tevatron, and the Large Hadron Collider at CERN. One of the most important objects is 
the parton distribution functions (PDFs). These functions describe the internal structure 
of hadrons in terms of the distribution of the longitudinal momentum fraction x carried by 
partons in the infinite momentum frame and the relevant QCD factorization has been well 
developed jlj. 

In recent years, hadronic physics community have extended the Feynman parton distri- 
butions to include the dependence on additional dimensions, in particular in the transverse 
directions perpendicular to the parenting hadron momentum direction. These extensions ap- 
pear in two different fashions: in the transverse coordinate space as the generalized parton 
distributions (GPDs) [2]; in the transverse momentum space as the transverse momentum 
dependent parton distributions (TMDs) [§]. A number of experimental facilities, such as the 
12GeV upgrade of Jefferson Lab, the Relativistic Heavy Ion Collider at the Brookhaven Na- 
tional Lab, and the planned Electron-Ion Collider, are trying to measure these distribution 
functions. These studies will lead us to the final answers to the important questions con- 
cerning the nucleon (nucleus) structure: the proton spin and parton saturation in nucleon 
(nucleus) at small x. 

Transverse momentum dependence in the parton distributions is also crucial to under- 
stand some novel hadronic physics phenomena in high energy scattering processes. This 
includes, for example, the single transverse spin asymmetries and small- x saturation 
phenomena (1,0). In the case for the small-x physics, the so-called /^-dependent parton dis- 
tributions contain resummation effects which come from multiple scattering associated with 
the nucleus target. Phenomeno logically, the k t - dependent gluon distribution (also called 
unintegrated gluon distribution) function has been applied to describe various high energy 
hadronic processes |9|. 

In general, the initial and final state interaction effects associated with the transverse 
momentum dependent parton distributions introduce additional QCD dynamics in these 
processes. For example, they lead to the non-universality for the transverse spin dependent 
TMD parton distributions (3^6, ID -15]. In ref. [l6j], we extended the universality discussions 
of TMD parton distributions to the small- x domain. In order to study the factorization 
issues relevant to small- x saturation physics, it is of advantage to focus on the two scale 
processes, such as semi-inclusive hadron production in deep inelastic scattering, Drell-Yan 
lepton pair production in hadronic reactions, and the di-jet correlations in these processes. 
In particular, we analyzed the small- a: transverse momentum dependent quark distributions 
probed in hadronic dijet-correlation in nucleon-nucleus collisions, as compared to that in the 
deep elastic lepton- nucleus (nucleon) scattering. Due to the nuclear enhancement, any soft 
gluon exchanges originated from the proton can be neglected. Thus, we only need to resum 
the soft gluon exchanges with the nucleus target in the large nuclear number limit. This 
procedure eventually helps us to obtain an effective kt factorization with the modified nuclear 
parton distributions in pA collisions. There have been interesting experimental results on 
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di-hadron correlation in deuteron-gold collisions at RHIC, where a strong back-to-back de- 
correlation was found in the forward rapidity region of the deuteron as compared to the 
narrow back-to-back peaks observed in the central rapidity region [171 ] . The purposes of 
this paper are to derive in details the QED results which we described in ref. (l6[ where 
we have shown the non-universality of TMD parton distributions, and then generalize to 
small- x models in QCD. 

We will start with a scalar-QED model calculation. There are a number of motivations 
for doing this. First, QED model is calculable, which makes it straightforward to resum 
all order initial and final state interaction effects in various processes. This allows us to 
rigorously discuss these effects, and shed light on the real QCD calculations. Second, there 
is similarity between the QED model calculations and the QCD saturation models. Quite 
a few results in the former framework can be directly generalized to the latter one. Third, 
the QED calculations are important in their own perspective. High order corrections in 
QED processes are interesting topics and have attracted intensive investigations since the 
QED was founded several decades ago. In particular, for the lepton pair production and 
photon radiation processes associated with large nucleus, hi gher order QED corrections have 
generated interests from both theory and experiment sides [l8j. Current running heavy ion 
collisions experiments are pursuing these studies at both RHIC and LHC facilities. The 
theoretical investigations shall provide further understanding of these processes. 

In the saturation domain of QCD, the McLerran-Venugopalan model [l9[ describes high 
density small- X partons in a relativistic large nucleus by treating the nucleus as a set of 
randomly distributed color sources p a (z~,z±) which generate soft classical gluon fields. 
Using the McLerran-Venugopalan model, we will demonstrate that the k t dependent quark 
distributions at small- x are not universal and the quark distributions involved in DIS and 
di-jet production processes are distinct. Furthermore, we find a simple formula to relate 
these two quark distributions through a convolution with the normalized unintegrated gluon 
distribution. 

The rest of the paper is organized as follows. In Sec. II, we construct a scalar QED 
model to investigate the universality property for the small- x parton distributions in various 
processes, where all gauge boson exchange contributions can be summed up, including all 
initial and final state interactions. In Sec. Ill, we demonstrate that the small- x TMD quark 
distributions of a large nucleus are not universal as well. In Sec. IV, we summarize our 
results. 



II. INITIAL/FINAL STATE INTERACTIONS IN QED MODELS 

In this section, we discuss the universality issue for the transverse momentum dependent 
parton distributions, by studying a scalar QED model. We will first introduce the QED 
model for our calculations, and discuss the universality issue of the TMD parton distribution 
functions in several processes. 



A. Model Description 

We follow Ref. [s| to construct the model for our calculations. This is a QED scalar 
model. It consists of heavy D and light (ft charged scalars with masses M and m, respectively, 
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(a) (b) (c) 

FIG. 1. Illustration of the final state interaction effects in DIS process (a), initial state interactions 
in Drell-Yan process (b), and initial/final state interactions in dijet-correlation process (c). In the 
dijet-correlation process, the gauge boson shall couple to any of the initial/final quarks. 

interacting with massive U(l) gauge fields with the mass A, 

Aqed = (P^X^i + tlW (1) 

via the covariant derivative = + igA^. In the above equation, we introduced two 
charged scalar particles, 0i and 02 with charges g\ and g>2, respectively. The purpose for 
this choice is to study the universality of the parton distribution in dijet-correlation. 

We will adopt this model to calculate the DIS, Drell-Yan, and dijet processes. In par- 
ticular, we will study the associated quark distribution functions in the small-x limit and 
investigate the universality issue. In Fig. 1, we plot the schematic diagrams for these pro- 
cesses in the scalar-QED model: (a) for DIS; (b) for Drell-Yan; (c) for dijet-correlation. In 
the DIS process (Fig. 1(a)), virtual photon scatters on the scalar quark from the nucleus 
target. In our example, the nucleus target has strong coupling with the Abelian gluon, 
g ^> 1. We need to resum all order gluon exchange contributions with the nucleus target, 
for which we have shown in the diagram. These interactions are referred as final state in- 
teractions. Similarly, we have multi-gluon interaction contribution between the incoming 
scalar quark with the nucleus target in the Drell-Yan process as shown in Fig. 1(b). For 
the dijet-correlation process in Fig. 1(c), there are both initial and final state interaction 
contributions. 

The quark distribution functions in the DIS and Drell-Yan processes in this scalar-QED 
model have been calculated in the literature. Let us recapture the main step in these cal- 
culations. In the high energy scattering process, we apply the power counting method to 
separate short distance physics from that from long distance. This effectively factorizes the 
cross section in terms of parton distributions. In the multi-gluon exchange contributions 
illustrated in Fig. 1, the dominant contribution in high energy limit comes from the parton 
distribution in nucleus target. For example, in the DIS process, the leading power contri- 
bution to the differential cross section can be factorized into the quark distribution from 
nucleus in Fig. 1(a). Of course, higher order corrections will be taken into account for the 
gluon radiation diagrams. 

The first step for this factorization is the eikonal approximation, which is valid in the 
leading power contribution, i.e., 1/Q 2 for DIS and Drell-Yan processes. Under this limit, 
the final state interaction contribution diagram can be simplified as the eikonal propagator, 
which can also be summarized as the gauge link contribution from the associated parton 
distribution definition in these two processes. For the dijet-correlation, since there is no 
simple definition, we will not seek the gauge link definition for that, although there has been 
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FIG. 2. Eikonal propagators for the initial/final state interactions in the scalar-QED model: (a) 
final state interaction on the scalar quark line; (b) final state interaction on the scalar antiquark 
line; (c) initial state interaction on scalar quark line. 



attempt to do that in the literature [111]. However, we emphasize in the leading power con- 
tribution, we shall be able to obtain the effective parton distributions in terms of the parton 
transverse momentum and longitudinal momentum fraction. We notice that higher order 
corrections will introduce large logarithms. To correctly resum these large logarithms, we 
need to pay special attention to the transverse momentum dependent parton distributions. 
At the current level of this paper, we do not need to worry about this additional effect. 

In the scalar-QED model, the eikonal approximation leads to the final state interaction 
effect in Fig. 2(a) as, 



which is the same as that for the fermionic propagator. Similarly, we will obtain the eikonal 
propagator for the final state interaction on the scalar antiquark as illustrated in Fig. 2(b), 
and the initial state interaction on the scalar quark in Fig. 2(c). The initial state interaction 
on the scalar antiquark will be opposite to that in Fig. 2(c). For the charge gi scalar quarks, 
we will have the same expressions by replacing gi with g 2 in these diagrams, respectively. 
In the following calculations, we will utilize these eikonal approximation for the relevant 
Feynman diagrams. 



B. Universality of Quark Distributions Between SIDIS and Drell-Yan Processes 

In this subsection, we review the known results 0,0, 2(j for the quark distributions in the 
above described model for the DIS and Drell-Yan lepton pair production processes in the 
small-x limit. In particular, we will perform the eikonal approximation on the final state 
interactions on the quark line. For the antiquark line with momentum p 2 , because its phase 
space is integrated out, we will keep the full kinematic dependence in these diagrams. The 
scalar quark (with charge g{) distribution in DIS process can be written as 

? DI W) = ^/ ^ / ^P + V^99,r^\k,M- q ,)\\ (3) 
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where the first three expansions of the amplitude are found as follows 
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(4) 
(5) 
(6) 



For convenience, we have defined the following integral, 

d 2 k l± d 2 k 2± 1 1 



d[l]d[2] 
d[l]d[2]d[3] 



(2vr) 4 k 2 1± + A 2 k 2 2± + A 2 
d 2 ki±d 2 k 2 ±d 2 k 3 ± 1 1 



(2tt) 2 5( 2 )(A; ± -A; 1± -A; 2± ) , 
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A; 2 ± + A 2 A; 2 ± + A 2 k 2 ± + A 2 



(2tt) "> 1± -r /v n, 2± 
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(7) 



We have also defined D (p±) = 2xP + p 2 + p\ L + m 2 , D\ = D(q±), D 2 = D(p 2 ±), Du = 
D(qi_ — ki±) and D 2i = D(p 2 j_ — ki±). To perform the resummation, we introduce the 
following Fourier transform, 



A(R ± ,r ± ) 



d 2 p 2 _i_ d 2 ki 



;A(k±,p 2± )e 



-ik ± -R ± —ip 2 ±-r± 



(2tt) 2 (2tt) 2 ' 

By applying this Fourier transform, we obtain, 

A$j s = -V(r ± )W(r ± ,R ± ) , 
Agl = +^V(r ± )W 2 (r ± ,R ± ) , 
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A£ls = +^V(r ± )W 3 (r ± ,R ± ) , 



where the functions V and W are defined as 

1 



V(r ± ) = —K (M r ± ) , 

Z7T 
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with Mq = 2xP + p 2 +m 2 , and Kq is the Bessel function. Clearly, we can see that the above 
expansion comes from the following exponential form, 



A DIS = iv( r± ) [1 - e- i99lW(r ^} 



(14) 



Substituting the above result into the quark distribution expression, we will obtain the TMD 
quark distribution in the DIS process, 
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In Sec. Ill, we will discuss how to translate this result into the fermion quark distribution 
function at small- x for this process in QCD and compare to the color-dipole/color glass 
condensate formalism. 



Similarly, one can calculate the quark distribution |20[ in the Drell-Yan lepton pair pro- 
duction process. The definition of the quark distribution has the same expression as in 
Eq. ([3]). The amplitudes of first three orders are found to be 
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(16) 
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It is easier to do the resummation when we perform the Fourier transform, for which we 
have 



A®=-V{r ± ) G(R±) - G(R± + f±) 



DY 

1(2) 
^DY 

(3) 



.4 
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where G(R±) = ^K (XR±). The all order resummation leads to 
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There is infrared divergence when the gluon mass goes to zero A — > 0. However, this infrared 
divergence cancels out when we calculate the quark distribution. Furthermore, taking the 
limit A 0, we also have —G(R± + r*j_) + G(R±) — >■ W(r±, R±), and we can find that the 
quark distribution is universal between DIS and Drell-Yan processes 



~DY 



-BIS 



(23) 



This universality is also guaranteed by the time-reversal invariance. The quark distributions 
in these two processes can be connected through time-reversal transformation, and time- 
reversal invariance can show that they are the same. However, this will not be the case 
when we compare the photon-jet correlation in Drell-Yan type process and dijet-correlation 
in DIS type process. We will carry out these calculations in future studies. 



C. Non-universality of Quark Distribution Between Dijet-Correlation and SIDIS/DY 
Processes 

To study the TMD scalar quark distribution in hadronic processes (e.g., pA collisions), 
we introduce the dijet-correlation process as illustrated in Fig. 1(c), 

p + A -> Jetl + Jet2 + X , (24) 
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FIG. 3. Lowest-order graphs for di-jet production in a hadron-hadron collision in the small-x limit. 
In these graphs, there is one soft gluon exchange with momentum k in addition to the hard gluon 
exchange. 



where the transverse momenta of these two jets are similar in size but opposite to each other 
in direction. For simplicity, we focus on the partonic channel qq' — > qq' for now in this paper 
and study the quark distribution. We expect that similar conclusions can be drawn as to 
other partonic channels and thus to the gluon distributions as well. In the ideal case, these 
two jets are produced back-to-back. However, the gluon radiation and intrinsic transverse 
momenta of the initial partons induce an imbalance between them. We are particularly 
interested in the kinematic region that the imbalance q± = P\± + P 2 ± is much smaller 
than the transverse momentum of the individual jet, namely, \q±\ <C |-Pi±| ~ |-P2±|, which 
also corresponds to the kinematics in the STAR measurements. Only in this region, can 
the intrinsic transverse momentum have significant effects. Since there are two incoming 
partons, both intrinsic transverse momenta can affect the imbalance between the two jets. 
For large nucleus and small-x, the dominant contribution should come from the intrinsic 
transverse momentum of the parton from the nucleus, which we label as q± in Fig. 1(a). In 
the following, we will focus on this contribution. 

Since we are interested in studying the final state interaction effects on the parton dis- 
tribution of the nucleus, for convenience, we choose the projectile as a single scalar quark 
with charge g 2 , which differs from the charge of the scalar quark from the target nucleus, 
g\. In addition, we assume that the Abelian gluon attaches to the target nucleus with an 
effective coupling g which is much larger than g2 or g\. All the partons in this calculation 
are set to be scalars with a mass m. The coupling g 2 being different from g\ is to show the 
dependence of the parton distribution on the initial/final state interactions associated with 
the incoming parton. If the dependence on g 2 remains for the nucleus parton distributions, 



they are not universal [131 . 114 



In Fig. EJ we plot the lowest-order graphs containing one soft gluon exchange with the 
momentum k. Following the discussions above, we keep the low transverse momentum 



approximation in terms of q±/Pi± (q±/p2±) by applying the power counting method [12 
Again, the important simplification is the eikonal approximation, which replaces the gluon 
attachment to the initial and final state partons with the eikonal propagator and vertex. Af- 
ter taking the leading order contributions, we find that the q± dependence of these diagrams 



can be cast into an effective quark distribution [121 ]. which takes the following form 
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FIG. 4. Example diagrams for two (a) and three (b) gluons exchanges, where the gluons can attach 
all charged particles in the upper part of the diagrams to the nucleus target. 



with p± = k± — q±. Here, the hard partonic part depending on the hard momentum 
scale Pi± has been separated from the above quark distribution in the differential cross 



section JjJ. This separation is only possible at the leading power contribution of q±/Pi 
The contributions from Fig. 2 can be written as, 
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99i 
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k 2 ± + A 2 2p~ 
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—k + + ie k + + ie 



(26) 



where Di follow the definitions introduced above. The first and second terms in the first 
square bracket correspond to Fig. [3] (a) and Fig. [3] (b), respectively. Fig. [3] (c) and Fig. [3] (d) 
yield the contributions as shown in the last term in Eq. f )26|) . It is not hard to see that these 
two contributions simply cancels since the sum is proportional to 2iri5(k + ) while k + ^ 0. 
This means that at the leading order in the coupling constant the dependence on g 2 drops out, 

which will however change at higher orders. Thus, one gets (k,p) = gg\ 



i 



k 2 , +A 2 



1 

Di 



1 

D 2 



At the next-to-leading order and the g 3 order, there are 20 and 120 graphs in total 
in covariant gauge, respectively. We show one of these graphs as an example in Fig. 
(a) and (b), respectively. Additional diagrams can be obtained by attaching the gluons 
to all incoming and outgoing scalar quarks. We organize our calculations according to 
these attachments. Again, the eikonal approximation discussed at the beginning of this 
sections will be utilized in the evaluations of these diagrams. For example, by applying 
this approximation, we replace the final state QED-gluon interactions with the scalar quark 
line with the eikonal propagators and vertices associated with g± coupling, which we label 
with rii hi Fig. |5j Similar approximations are made for all other diagrams, with n 2 and 
representing final and initial state QED-gluon interaction with the scalar quark line with 
charge g 2 - We summarize those diagrams in Fig. (J5]), where every graph except for (e) and 
(f) represents 3 graphs and total of 20 graphs contribute. 

Among those 20 graphs, only a few sets of graphs give non-vanishing contributions. To 
better explain the calculation, let us first evaluate Fig. (JSJ) (a) and (b) as an example. 
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FIG. 5. Two-gluon-exchange contributions to the quark distribution, z = 1,2,3 for each m in each 



graph with n\ 
(h). 
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_ k 9 + +it and n 3 = ^Bpfe- * ^ 3 for n i,j is implied in graphs (g) and 



Fig. ©(a) and Fig. ©(b) with rii = rt\ give the following contribution 
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After integrating over dk^ and dfc^~, one finds that the subtotal contribution reads 

111 
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It is straightforward to find that the contributions from n 2 and cancels and g 2 contribution 
vanishes due to the same reason we explained under Eq. (1261) . 

As to the graphs illustrated in Fig. (j5j) (c) and Fig. ©(d), we find that the contribution 
from ns vanishes due to vanishing contour integral of dk£. Thus, the subtotal contribution 
of Fig. ©(c) and Fig. ©(d) is 
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with k = ki + k 2 . After integrating over d/c^~ and dk 2 , one finds 
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Similarly, we obtain 
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The total contributions are 
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where J d[l]d[2] stands for the same definition as in Eq. (j7|), D Xi = D(q A 
D(p± — kij_). Clearly, the second order result shows a dependence on g 2 . However, in 
the amplitude squared calculation for the quark distribution Eq. (|25|) . the g 2 dependence 
from A^'A^ 2 ' is canceled out by its complex conjugate because A^ is real while A^ 
is purely imaginary. The leading order contribution to g 2 comes from \A^ (k,p)\ and 
A^A^* + A^A^*. Therefore, to see the residue dependence on g 2 , we need to carry out 
the calculation of the amplitude up to order g 3 . 

At the g 3 order, there are 120 diagrams in total with three soft gluon-exchange (see e.g., 
Fig. HI (b)), including all possible permutations of the attachments of these three gluons to 
the target nucleus. Let us take Fig. H] (b) as an example together with the other 5 crossing 
diagrams. The corresponding contribution is 
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In reaching above result, we have used the following two integrals: 
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Summing up all these graphs, we obtain the three gluon exchange amplitude, 
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where J d[l]d[2]d[3] follows the same definition as in Eq. ((?]). Again, we see the dependence 
on g 2 in the second and third terms. An important cross check of these results is that, if we 
set g 2 = ~9i, there is effectively no charge flow in the final state, and the quark distribution 
is identical to that in the Drell-Yan process in the same model. Applying g 2 = —gi, we can 
easily see that indeed we reproduce those calculated in Ref. |20J. Also, by setting g 2 = 0, 
we can recover the DIS amplitudes. 

With the amplitude calculated up to O (g 3 ), we are able to check the dependence on g 2 
for the parton distribution in Eq. (125]) . Substituting the above amplitudes into Eq. (|25|) . we 
find that the g 2 dependence still remains up to order g A . If we drop all g 2 terms in these 
results, we obtain the quark distribution in DIS in the same model |6|, |8j . This clearly shows 
that the TMD quark distribution q(x, q±) is not universal. 

This non-universality is better illustrated when we sum up all order multi-gluon exchange 
contributions. To do that, we introduce the following Fourier transform jsf, 



A(R,r) = J 



d 2 k± d 2 p± 

(27T) 2 (27T) 2 



e~ ik± - R± ~ ip± - r± A(k,p) . 



(38) 



From the Fourier transforms of A ( - 1 ' 2 ' 3 \k,p), we can easily see that they follow the expansion 
of an exponential form, 



A** ** (R, r) = Y, A{n) (R, r ) = iV (r ± ) {l - e 



iggi[G(R x +r ± )-G(R ± )] 



} 



-i992G(R ± ) 



(39) 



where G(R±) = K (XR±) /2vr and V(r±) = K (Mr ± ) /2vr with M 2 = 2xP + p- + m 2 . In 
the above result, the ^-dependence seems to only appear as a phase which may not lead to 
a physics consequence. However, because the transverse momentum qj_ is conjugate to the 
coordinate variable difference R±-r±, this phase will lead to a non-universal contribution to 
the quark distribution as defined in Eq. (I25j) . Here we can identify the factor e~ t99iG ^ R± ^ 
as a Wilson line which essentially resums soft interactions between quarks and the target 
nucleus. 

Therefore, the all order result reads, 



xP+ 2 



J dp-p- J d 2 R ± d 2 R' ± d 2 r ± d 2 r' L 5 (2) (R ± +r±- R' ± - 

Xe -iqx{r x -r' ± ) e -ig g2 (G(Rs_)-G(R' ± )) v ^ 

x jl - e i99i[G{R ± +r ± )-G(R ± )}^ jj _ ^iggy [g{& x -*> x )-Q(b! x )\ 



(40) 



This TMD quark distribution is clearly different from that calculated in DIS in the same 
model 0,0]. In other words, TMD quark distributions are not universal. It is interest- 
ing to notice that the g 2 dependence disappears after the integration over the transverse 
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momentum. This is consistent with the universality for the integrated parton distribu- 
tions @, 0, E3- In spite of the non-universality, we expect that one can still reach an 



effective TMD factorization formula for pA collisions by absorbing all the violation effects 
into the parton distributions as in Eq. fl40p . 

It has been argued that the light-cone gauge may simplify the factorization property for 
the hard scattering processes. For example, if we choose the advanced boundary condition 
for the gauge potential in light-cone gauge, the wave function of hadrons contain the final 



state interaction effects [6|, |21|. However, as we showed in the above calculations, this does 
not help to resolve the (^-dependence in the quark distribution in the dijet correlation in 
hadronic process due to the presence of both initial and final state interactions. In other 
words, the quark distribution from the nucleus target has to contain the interaction with 
the incoming (outgoing) quark with charge g2, which can not be solely included into the 
wave function of the nucleus target. 



III. EXTENSION TO THE QUARK DISTRIBUTIONS AT SMALL-x IN CGC 



In this section, we extend the previous calculation to the quark distribution functions 
for a large nucleus at small- x in QCD, by calculating the similar resummation effects due 
to initial and final state interactions. As an example, we follow the McLerran-Venugopalan 
(MV) Model [lH . The MV model describes high density gluon distribution in a relativistic 
large nucleus by solving the classical Yang-Mills equation. An effective theory, called the 
Color Glass Condensate (CGC), is developed to study the high density physics in QCD in a 
systematic manner It is equivalent to the saturation picture which is based on the color 



dipole model [22j, |23j in terms of quantitatively describing the parton saturation at small-x. 



In our following discussion, we will focus on the CGC formalism since it is close to what we 
have used above in the QED model. 

Following the CGC formalism, the target nucleus is treated as a collection of randomly 
distributed color sources p a (z~, z±) which generate soft classical gluon fields. Similar to 
what happens in above QED model, where we obtain the Wilson line e~ l99i ^ G ^ x± ^, the soft 
interactions (i.e., the soft gluon exchanges) between the nucleus and a relativistic quark can 
also be resummed into a Wilson line which reads 



U (x±) = Texp 



~W9i 



J dz d 2 z±G (x± — z±) p a (z , z±) t a 



(41) 



where t a is the SU(3) color matrix in the fundamental representation and the two- 
dimensional propagator G (x± — z±) is the same as the one we used in the scalar QED 
model. The differences come from the fact that the target nucleus is no longer treated as 
point particles and quarks now carry colors. It is straightforward to see that we recover 
the scalar-QED model result (U (x±) =>• e~ iggi ^ G( - x± ^ if we apply the point particle approx- 
imation: p a (z~,z±)t a S (z~~) 5^ (z±). In addition, the ensemble average over the color 
sources should be performed. Since one assumes the color sources are randomly distributed, 
a Gaussian distribution W [p] is always used in the average. In the McLerran-Venugopalan 
Model, the Gaussian distribution is defined as follows [s| 



W [p] = exp 



dz d 2 z± 



Pa{z ,Z±)p a (z ,Z ± ) 

2f (z~) 



(42) 
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The variance p 2 (z~) of the charge distribution represents the density of color sources per 
unit volume. It is assumed that the average in CGC is performed in a functional integration 
over p accompanied by W [p] . It then follows that 

(p a (x~,x±) p b (y~,yi_))p = J VpW[p]p a (x~,x ± ) p b (y~,yA_) (43) 

= S ab S (x- - y-) 6® (x ± - y ± ) p 2 (x-) . (44) 

This turns out to be useful in the following derivations for the TMD parton distributions in 
CGC. 



A. DIS and Drell-Yan Processes 

As we discuss in Sec. II, the TMD quark distribution in DIS in the scalar QED model can 
be written as 

q B1S (x,q±) = ^ / dp-p- [ d 2 R ± d 2 R' ± d 2 r ± d 2 r' ± 

x 

X 



<5 (2) {R ± + r ± - R' ± - r' ± ) e tq ^- R '^V{r ± ) V {r' ± ) 

|l _ e iggi[G(R ± +r ± )-G(R ± )]^ | j _ e - lggi [G(R' ± +r' ± )-G(R' ± )}^ _ ^ 



To derive the above result, we have assumed the target hadron (nucleus) to be a point 
particle. In order to calculate the quark distribution in CGC, we need to relax the point 
particle approximation. Following the above discussions, we first assume that the target 
hadron has a color charge distribution p a (z~, z±) and perform a replacement e~ t991 ^ G ^ x± " =>• 
U{x x ). 

The second step is to average over the color sources p a (z~,z±_), which appears in the 
exponents of the Wilson lines U (x i ), with the Gaussian distribution W [p\. Following this 



procedure, one finds (see Refs. |24j-|26[|) 

e i 99l [G(R ± +r ± )-G(R ± )] ^ Tr ^ t ( R± + r± ) JJ ( R± )j 



= iV c exp |-/i 2 J d 2 z ± [G {R ± + r± -z ± )-G {R x - z ± )} 2 ^ 

= N c exp[-p 2 s J dV[G(rj_ + zJ-G(zjJ] 2 J 

~ N c ex P {-Q 2 s r 2 j4}, (46) 

where N c = 3 is the number of colors, the saturation scale Q s is defined as Q 2 = ^ In -^-^ 

with p 2 = ^-^t a t a jdx~p 2 {z~). In the evaluation of the above two-point functions 
(U* (R± + r±) U (R±)) , we have assumed that the nucleus size is so large that we can 
shift R± in the transverse integration. The saturation momentum naturally arises as a 
result of multiple scatterings between the hard parton and color charges inside the nucleus. 

The next step is to use fermionic quark splitting kernel instead of the scalar quark splitting 
kernel. Thus we replace V (r±) in Eq. (|4"5l) by k-2Ki (Mr±) where the factor of 2 comes 
from the fact that fermionic quark has two different helicities. It is straightforward to derive 
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this fermionic quark splitting kernel as in Ref. 27|. The rest of the calculation will remain 
the same since the eikonal propagator for a fermionic quark is the same as the one for a 
scalar quark as in Eq. (j2j). After changing the integral variable to y = 2xP + p~, we can cast 
the quark distribution into 



xq ms (x, q ± ) 



32tt 6 



dyd z R±d z rj_d 2 r' ± e 



Vr ± K (y/yr±) ■ X7 r > K (y/yr' ± 



x 



1 + exp 



Ql (r ± - r' ± 



exp 



qy 



4 



exp 



2^/2 



• (47) 



The virtuality of the virtual photon Q 2 = 2xP + P~ is taken to be much larger than Q 2 and 
q\. Therefore, one can approximately integrate y from to +oo. The dominant contribution 
comes from the region where y is close to 0. It is hard to evaluate above integrals analytically. 
Nevertheless, we can study the quark distribution in the large and small q\ limit, which give 



dxg DIS (x, q±) 



d 2 R ± 
dxg DIS (x, q±) 



d 2 R ] 



g 2 t »Qi 



q 2 <Q2 



N c Q 2 
12vr 4 q 2 ± 



4vr 4 



(48) 



These results agree with those derived in the saturation model for the quark distribution of 
a large nucleus in DIS (see e.g., Eqs. (27-29) of Ref. @fl Furthermore, we can transform 
the above results to the momentum space and define the normalized unintegrated gluon 
distribution F(k±, Q s 



as 



F{k ± ,Q s ) 



d 2 r±_ 

WY 2 



-ik i 



Tr(U(R ± )W(R ± + r ± )) p 1 



exp 



Q 2 



IT, ■ (49) 



In arriving at the Gaussian form of F(k±, Q s ) in Eq. (149 p . we have neglected the logarithmic 
dependence of r\ in the saturation momentum Q s . Thus, one can write the quark distri- 
bution as a convolution of the unintegrated gluon distribution and the splitting kernel in 
momentum space, 



xq ms (x, q ± ) 



4tt 4 

Nc_ 
4vr 4 



d'R ± d z k ± F(q ± - k ± , Q s ) / dy 



q± 



k 



d 2 R ± d 2 k ± F{k ± ,Q & 



Q± 



q± + y 
(q± - k± 



k 2 ± + y 



In 



ql 



(50) 



which is consistent with the results obtained in Ref. 28, 29||. The unintegrated gluon distri- 
bution F(k±, Q s ) is usually defined through the scattering amplitude of a dipole with size 



r± on the target nucleus 30, 3l[. This dipole scattering amplitude is also equivalent to the 
expectation value of a Wilson loop with width r± and infinite length as we used above. 



1 We notice that there is a factor of 1/2 difference between our results and those obtained in Ref. 27|. This 
difference comes from the fact that the quark distribution calculated in Ref. 27[ is in fact the total quark 
distribution which includes anti-quark distribution as well. 
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For comparison, we can also calculate the quark distribution involved in the Drell-Yan 
process. Again, we start with scalar-QED model result |20j |. 



g DY (x,g ± ) 



x 



p+2 



dp-p- I d' z R±d 2 R x d 2 r±d z r 



x 



<5 (2) (R ± + r ± -R' x - r' ± ) e iq± ' 



, R±—R', 



>V(r ± )V(r' ± 



x < e 



ig 9l G(R ± +r ± ) _ iggiG(R ± )\ f -iggiG(Rf ± +r' ± ) _ -ig 9l G(R' ± ) 



(51) 



Following the same procedures, we find that the quark distribution in the Drell-Yan process 
in the CGC formalism reads as 



xq DY (x, q ± ) 



Nr. 



Jdyj d 2 R ± d 2 r ± d 2 r' ± e-^i r ^) Vr± K (Jyr ± ) • Vr > ± K (y/yr' ± ) 
-Ql (r± ~ Af /A] - exp [-QyjA] ~ exp [-Q 2 s r' 2 jA\ } . (52) 



32vr 

x <^ 1 + exp 



The quark distribution in the Drell-Yan process is the same as that in DIS, which is consistent 
with the conclusion in the scalar QED model, as the QCD factorization predicts. 



B. Dijet production in pA Collisions 

Finally, let us consider the TMD quark distribution for a large nucleus involved in the 
di-jet production, again, taking the qq' — > qq' channel as an example. In the scalar-QED 
model, the TMD quark distribution in this process is shown in Eq. (HOl) . In order to extend 
to the real QCD calculation, we will assume that the color charge for quark q' is the same 
as the quark q in the sense of the average over the large nucleus. This means that we will 
set $2 — 9\ i n the scalar-QED result El- Furthermore, we find that the quark distribution 
here will naturally involve four-point function. For example, expanding the phase factor in 
Eq. fj4T)]) will depend on the four-point function in the CGC formalism, 



(53) 



iggi(G(R x )-G(R' x )) l\ _ iggi[G(R x +r ± )-G(R x )]\ _ -iggi[G(R' x +r' ± )-G(R' x )] 



U (R ± ) W (R' ± ) + U (R ± ) W (R' ± ) U (R ± ) W (R' ± ) 
-U (R ± ) W (R' ± ) U (R ± ) W (R ± +r ± )-U (R ± ) W (R' ± ) U (R' x + r' ± ) W (R' ± 



The appearance of the four point functions signals the difference between the di-jet produc- 
tion process and DIS, whereas the latter only involves two point functions. This indicates 
that parton distributions directly extracted from DIS are not sufficient to compute and 
describe the dijet production processes. 

We can further simplify the above result by taking the large N c limit for the four point 
functions [32| 

(U (R ± ) f/ f (R' ± ) U (R ± ) C/t (R ± + r ± )) p ~ exp j-^ [(r ± - r' ± ) 2 + r 2 ± ] } . (54) 



2 This can be checked against the lowest nontrivial order perturbation expansion of the multi-gluon exchange 
contributions in the large N c limit. As we mentioned before, setting 92 = ~9i will lead to the quark 
distribution in the Drell-Yan process. 
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With this reduction, we arrive at the following quark distribution in the large N c limit, 



xq BJ (x,q±) 



327T 6 



dyd R±d r±d r ± e 



Vri K (y/yr ± ) ■ Vr' K (y/yr'±) 



exp 



r x -r _[ 



— exp 



;((rx-rl) 2 +rj) 



+ exp 
— exp 











to 




Ql((r ± -r' ± ] 


2 +r'l) 



(55) 



which then yields 



dxq (x,q±) 



d 2 R 
dxg DJ 



x,q± 



d?R ] 



A»Ql 12714 & 



0.44 



4vr 4 



(56) 



It is straightforward to see that the quark distributions in DIS and di-hadron production 
have the same perturbative tails while they differ in the small q\ limit. As shown in Fig. |6l 
the quark distribution is about twice broader than the one in DIS while its peak is about half 
of the peak of the DIS distribution. However, it is easy to check analytically and numerically 
that the integrated quark distributions are universal for these processes. 

In the momentum space, we find that the quark distribution in di-jet production can be 
written as follows: 



x~q m 



x,q±) 



Nc_ 
4vr 4 



X 



J d 2 R ± J d 2 l ± F(q ± -l ± ,Q s ) 

l± ■ (Z ± - kj 



d 2 k ± F(k ±} Q s 



In 



il - (i± - k L y (i ± - k ± ) 2 



which implies that 



xq m (x, q ± ) = j d 2 l ± xq DlS (x, l ± ) F(q ± - l ± , Q s ). 



(57) 



(58) 



This is an interesting new result. It relates the two apparently different quark distribu- 
tions through a k t convolution with the unintegrated gluon distribution F(q± — l±,Q s ). It 
is easy to see that both quark distribution reduce to the same form after integration over 
q± since F(q±_ — l±,Q s ) is normalized to 1. In addition, Eq. (j5SJ explains the broadening of 
the di-jet quark distribution as shown in Fig. O This formula has a natural physical inter- 
pretation. This convolution arises as a result of the extra initial and final state interactions 
in the di-jet production process. 



IV. SUMMARY AND DISCUSSIONS 

In this paper, we have studied the initial and final state interaction effects in the small-x 
parton distributions. As an example, we discussed the quark distributions in the semi- 
inclusive deep inelastic scattering, Drell-Yan lepton pair production and dijet-correlation 
in pA collisions. We calculated these distributions first in a scalar-QED model and then 
extended to the CGC formalism in QCD. 
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2 4 6 8 



q t 2/c v 



FIG. 6. Comparison of quark distributions a s functions of |p in DIS (or Drell-Yan) 

and di-hadron production. The solid curve stands for the quark distribution in DIS and Drell-Yan 
process, and the dash curve represents the distribution involved in di-hadron production. 



We have shown the non- universality for the small- x parton distributions in dijet correla- 
tion in the scalar QED model calculations, as compared to the quark distributions in DIS 
and Drell-Yan processes. For the particular partonic channel qq' — > qq', we find that the 
net effects are summarized into a phase which leads to a non-vanishing contribution to the 
quark distribution and breaks the universality. 

We have also calculated the TMD quark distribution involved in dijet production in the 
saturation models 27|| • We reached the conclusion that TMD quark distributions are not 



universal in the color-dipole or color glass condensate formalism, by showing that the quark 
distributions involved in DIS and di-jet production processes are distinct as shown in Fig. |6j 
In addition, we found a simple formula which relates these two different quark distributions 
through a convolution with a normalized gluon distribution. 

It is interesting to note that we can also compare to the quark distribution discussed 



in Refs. 11 , where the initial and final state interaction effects are summed into an ef- 



fective gauge link associated with the quark distribution. For example, for the partonic 
channel qq' — > qq', the quark distribution requires the gauge link as G = jf^—^ — -WW — 

j-W^WH which is different from that in the semi-inclusive DIS process with G = 
In the large N c limit, the additional gauge link structure would contribute a factor which is 
similar to the unintegrated gluon distribution in the CGC formalism as we have shown in 
Sec.III. 

Despite the non-universality, we expect that there exists a generalized TMD factorization 
for the di-jet production in pA collisions in the large A limit. Thanks to the nuclear en- 
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hancement, which allows us to neglect any soft gluon exchanges originated from the proton, 
we can resum all the anomalous terms which breaks the k t factorization and put them into 
the parton distributions of the target nucleus. This procedure leads to an effective k t fac- 
torization with non-universal nuclear parton distributions in pA collisions. Also we would 
like to emphasize that the Wilson line U(x±), which provides the underlying fundamental 
description of the interactions between partons and dense hadronic matter, is still universal. 

Of course in the high energy limit, the forward di-jet production in pA collisions is 
dominated by the qg — > qg channel due to high gluon density in the target nucleus. This 
calculation is much more complicated than the qq' — » qq' process which we considered in 
this paper. The complexity comes from the fact that there are many more channels involved 



in the qg — > qg process. There have been some theoretical calculations [33l. l34j | in CGC. 



However, the non-universality issue has not yet been taken into consideration. We will 
address this problem together with the photon-jet productions in pA collisions in a future 
publication [351 ]. 

The non-universality for the TMD parton distributions at small- £ clearly imposes a chal- 
lenge in explaining the dijet-correlation data in dA collisions at RHIC with the parton 
distributions extracted directly from the DIS data. The non-universality, on the other hand, 
provides an opportunity to study QCD dynamics associated with the initial and final state 
interaction effects, which are calculable at small- x (high gluon density limit) according to 
our results. More phenomenological discussion will be provided in ref. (35| . 



V. ACKNOWLEDGMENT 



We thank Tony Baltz, Stan Brodsky, Fabio Dominguez, Paul Hoyer, Cyrille Marquet, 
Larry McLarran, Al Mueller, Jianwei Qiu and Raju Venugopolan for interesting discus- 
sions. This work was supported in part by the U.S. Department of Energy under contracts 
DE-AC02-05CH11231. We are grateful to RIKEN, Brookhaven National Laboratory and 
the U.S. Department of Energy (contract number DE-AC02-98CH10886) for providing the 
facilities essential for the completion of this work. 



[1] J. C. Collins, D. E. Soper and G. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988). 
[2] K. Goeke , M.V. Polyakov, and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 (2001); 

M. Diehl, Phys. Rep. 388, 41 (2003); X. Ji, Ann. Rev. Nucl. Part. Sci. 54, 413 (2004); 

A.V. Belitsky and A.V. Radyushkin, Phys. Rep. 418, 1 (2005); S. Boffi and B. Pasquini, Riv. 

Nuovo Cim. 30, 387 (2007). 
[3] X. Ji, Phys. Rev. Lett. 91, 062001 (2003); A. V. Belitsky, X. Ji and F. Yuan, Phys. Rev. D 

69, 074014 (2004). 

[4] S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. B 530, 99 (2002); Nucl. Phys. B 642, 
344 (2002). 

[5] J. C. Collins, Phys. Lett. B 536, 43 (2002). 

[6] X. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002); A. V. Belitsky, X. Ji and F. Yuan, Nucl. 

Phys. B 656, 165 (2003). 
[7] D. Boer, P. J. Mulders and F. Pijlman, Nucl. Phys. B 667, 201 (2003). 



19 



S. J. Brodsky, P. Hoyer, N. Marchal, S. Peigne and F. Sannino, Phys. Rev. D 65, 114025 
(2002). 

E. Iancu, A. Leonidov and L. McLerran, arXiv:hep-ph/0202270; E. Iancu and R. Venugopalan, 
|arXiv:hep-ph/0303204[ J. Jalilian-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56, 
104 (2006); F. Gelis, E. Iancu, J. Jalilian-Marian and R. Venugopalan. larXiv:1002.0333l [hep- 
ph]; and references therein. 

D. Boer and W. Vogelsang, Phys. Rev. D 69, 094025 (2004). 

C. J. Bomhof, P. J. Mulders and F. Pijlman, Phys. Lett. B 596, 277 (2004); Eur. Phys. J. C 
47, 147 (2006); A. Bacchetta, C. J. Bomhof, P. J. Mulders and F. Pijlman, Phys. Rev. D 72, 
034030 (2005); C. J. Bomhof and P. J. Mulders, JHEP 0702, 029 (2007). 
J. W. Qiu, W. Vogelsang and F. Yuan, Phys. Lett. B 650, 373 (2007); Phys. Rev. D 76, 
074029 (2007). 

J. Collins and J. W. Qiu, Phys. Rev. D 75, 114014 (2007). 

W. Vogelsang and F. Yuan, Phys. Rev. D 76, 094013 (2007); J. Collins, arXiv:0708.4410l 
[hep-ph] . 

T. C. Rogers and P. J. Mulders, Phys. Rev. D 81, 094006 (2010) |arXiv: 100 1.29771 [hep-ph]]. 
B. W. Xiao and F. Yuan, Phys. Rev. Lett. 105, 062001 (2010) |arXiv: 1003.04821 [hep-ph]]. 

E. Braidot, for the STAR Collaboration, larXiv: 1008.39891 [nucl-ex]; B. Meredith, for the 
PHENIX Collaboration, to appear; 

see, for example, K. Hencken et al, Phys. Rept. 458, 1 (2008), and references therein. 

L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994) [arXiv:hep-ph/9309289] ; 

Phys. Rev. D 49, 3352 (1994) [arXiv:hep-ph/9311205] . 

S. Peigne, Phys. Rev. D 66, 114011 (2002). 

S. J. Brodsky, B. Pasquini, B. W. Xiao and F. Yuan, Phys. Lett. B 687, 327 (2010) 

|arXiv:1001.lT63l [hep-ph]]. 

A. H. Mueller, Nucl. Phys. B 415, 373 (1994). 

A. H. Mueller, |arXiv:hep-ph/9911289| 



F. Gelis and A. Peshier, Nucl. Phys. A 697, 879 (2002) |arXiv:hep-ph/0107142 1. 

J. P. Blaizot, F. Gelis and R. Venugopalan, Nucl. Phys. A 743, 57 (2004) 

|arXiv:hep-ph/0402257| . 

K. Fukushima and Y. Hidaka, JHEP 0706, 040 (2007) [arXiv:Q704.2806l [hep-ph]]. 
A. H. Mueller, Nucl. Phys. B 558, 285 (1999). 

L. D. McLerran and R. Venugopalan, Phys. Rev. D 59, 094002 (1999); R. Venugopalan, Acta 
Phys. Polon. B 30, 3731 (1999). 

C. Marquet, B. W. Xiao and F. Yuan, Phys. Lett. B 682, 207 (2009). [arXiv: 0906 .14541 [hep- 
ph]]. 



M. Braun, Eur. Phys. J. C 16, 337 (2000) arX iv:hep-p h/0001268|; Phys. Lett. B 483, 105 
(2000) [arXiv:hep-ph/0003003] . 

D. Kharzeev, Y. V. Kovchegov and K. Tuchin, Phys. Rev. D 68, 094013 (2003) 
|arXiv:hep-ph/0307037| . 

F. Dominguez, C. Marquet and B. Wu, Nucl. Phys. A 823, 99 (2009) |arXiv:0812. 38781 [nucl- 

th]]. 

C. Marquet, Nucl. Phys. A 796, 41 (2007); J. L. Albacete and C. Marquet. larXiv: 1005.4 065 
[hep-ph] . 

K. Tuchin, Nucl. Phys. A 846, 83 (2010) [arXiv:0912.5479l [hep-ph]]. 
F. Dominguez, B. W. Xiao and F. Yuan, to be published. 



20 



